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Abstract 

We estimate CMB polarization and temperature power spectra us- 
ing WMAP 5-year foreground contaminated maps. The power spec- 
trum is estimated by using a model independent method, which does 
not utilize directly the diffuse foreground templates nor the detec- 
tor noise model. The method essentially consists of two steps, (i) 
removal of diffuse foregrounds contamination by making linear combi- 
nation of individual maps in harmonic space and (ii) cross-correlation 
of foreground cleaned maps to minimize detector noise bias. For tem- 
perature power spectrum we also estimate and subtract residual unre- 
solved point source contamination in the cross-power spectrum using 
the point source model provided by the WMAP science team. Our 



1 



TT, TE and EE power spectra are in good agreement with the pub- 
lished results of the WMAP science team. The error bars on the 
polarization power spectra, however, turn out to be smaller in com- 
parison to what is obtained by the WMAP science team. We perform 
detailed numerical simulations to test for bias in our procedure. We 
find that the bias is small in all cases. A negative bias at low / in 
TT power spectrum has been pointed in an earlier publication. We 
find that the bias corrected quadrupole power (1(1 + 1)Ci/2tt) is 532 
fj-K 2 , approximately 2.5 times the estimate (213.4 fJ>K 2 ) made by the 
WMAP team. 



1 Introduction 



The anisotropics of the Cosmic Microwave Background (CMB) radiation are 
the most important evidence behind the tiny fluctuations tha t are gener- 
ated by the inflati o nary p aradigm of the Big-B ang cosmology (jStarobinsky 
( 1l982h : lGuth fc Pil (jl982h : lBardeen et all (jl988h ). One can d etermine cosmo 



l ogical par ameters precisely by measuring these anisotr opies (IJungman et al. 



( 11996bl Jal); iBond et al.l (119971 ); IZaldarriaga et al.l (119971 )). These anisotropies 
possess a certain degree of linear polarization due to the quadrupolar tem- 
perature pattern seen by the mov ing electrons in the primordial plasma (iRees 
(1l968h : lBasko k Polnarevl (Il98fjf )l Recently, the WMAP satellite has mapped 
the total intensity and polarization of the CMB anisotropies over the full sky 
in its 5 frequency b ands from 23 GHz to 9 4 GHZ with unprece d ented resolu- 
tion a n d sensitivity dBennett et al. (l2003al Jbl): [Page et all (!2007h : lKogut et al 



( 2007 ); Hinshaw et al. ( 20081 )). The polarization power spectrum acts as a 
complement to the temperature power spectrum. It leads to better con- 
straints on the cosmological parameters and is also useful to break degenera- 
cies among certain co smological par ameters, e.g. epoch of reionization and 
scalar to tensor ratio (iKmneyl (119981 )). Furthermore it h as been argued that 
CMB polarization may serve as a direct probe of inflation (jSpergel fc Zaldarriaga 
( 119971 )) , can test i f the p arity symmetry is pre served on the cosmological 



scales flLue et all (Il999h : iKomatsu et~ai1 (12008^ . can provide informatio n 



about the epoch when the first stars begin to form (ICrittenden et al.l (119931 ); 



Ng fc Ng ( 19951 )) an d provide a measure o f the g r avity waves that are gen - 



erated by inflation (lHarari &; Zaldarriagal (119931 ); ICrittenden et al.l (119951 ); 
Kamionkowski fc Kosowsky ( 19981 )). The WMAP team have produced their 



temperature and polarization power spectrum based upon the foreground 
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cleaned maps which are ob tained using p rior models of the s ynchr otron, dust 
and free-free components ( Kogut et al. ( 2007 ); Page et al. ( 20071 )). Though 
this method allows one to use all the available information about the fore- 
ground components it is also a very important scientific task to perform 
an independent analysis of the data by techniques which do not rely upon 
explicit foreground modeling. 

A multipole based approach for foreground removal was first proposed 
Tegmark fc Efstathioul ( 19961 ) and was implemented on the WMAP data 



by 
by. 



Tegmark et al.l (120031 ). Later, ISaha et al.l (120061 . 120081 ); ISouradeep et al. 



( 120061 ) extended this method to extract the temperature anisotropy power 
spectrum of the Cosmic Microwave Background (CMB) radiation from the 
raw WMAP data. The power spectrum is obtained by forming several cleaned 
maps using subsets of the available maps and thereafter cross-correlating the 
resulting maps. This Internal Power Spectrum Estimation (IPSE) method 
utilizes CMB data as the only input without making any explicit modeling of 
the diffuse galactic foreground components or detector noise bias. The fore- 
ground components are removed using the fact that in thermodynamic tern- 



(Mather et al. 


( 


1994): 


Fixsen et al. 



( 119961 ) ) . while the foreground components are frequency dependent. The de- 
tector noise bias is removed by cross-correlating different foreground cleaned 
maps obtained by using independent detector subsets. This substantially re- 
moves the noise bias since, to a good appr oximation, WMAP detecto r noise 
is uncorrelated for t wo different detectors (jJarosik et al.l (120031 . 120071 )). The 
final power spectra (ISaha et all (120061 . \2Q0m obtained by IPSE agrees well 
with the results published by the WMAP science team. Thus the method 
serves as an independent technique to verify the main power spectrum re- 
sult obtained by the WMAP science team starting from the stage of diffuse 
foreground components removal. The method has several advantages. First, 
the foreground components removal method is entirely independent of the 
foreground template models. Therefore the foreground cleaned maps are not 
susceptible to systematic errors that might arise in template based methods 
due to incorrect template modeling. Second, the cleaned power spectrum can 
be studie d analytical l y in th e special case of full sky one iteration foreground 
cleaning (ISaha et al.l (120081 )). This allows us to quantify and understand the 
statistical properties of the residuals in the cleaned power spectrum. This 
may be very useful in the case of noisy data or when the total number of 
available frequency bands are less than the total number of independent 
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parameters required for satisfactory modeling of all dominant underlying 
components. A detailed analytical s t udy o f the bias in the cleaned power 



spectrum is presented in IS aha et al.l (120081 ). Third, it is possible to obtain 
a model independent estimate of the map and powe r spectrum o f the t otal 
foreground emission at each of the frequency bands ( Ghosh et al. ( 20091 )). 

In the present paper we extract CMB polarization EE, TE as well as the 
TT power spectra using the WMAP 5-year foreground and detector noise 
contaminated maps as input to IPSE - our combined foreground removal 
and power spectrum estimation procedure. Since the CMB polarization sig- 
nal is weak the polarized maps published by the WMAP science team are 
dominated by the foreground components and detector noise. This seriously 
limits the accuracy with which the polarized CMB power spectrum can be 
extracted. However, since our method does not use any template model to re- 
move foreground components we argue that our power spectrum is free from 
systematic effects that might arise due to incorrect modeling of polarized 
(and temperature) foreground templates. 

The error bars as well as the bias in the extracted polarization power spec- 
trum are estimated by numerical Monte Carlo simulations. Here we make use 
of explicit foreground and detector noise models. In the case of temperature 
power spectrum a similar analysi s reveals the presence of a negative bias at 
low multipoles (ISaha et al. The bias corrected temperature spectra 

explains almost all of the low power observed in case of quadrupole. In this 
case it is also possible to analytically obtain an estimate of the bias in some 
special cases. 

Alternate approaches of CMB power spectrum estimation have been stud- 
ied by several authors, e.g., using foregrou n d clea n ed maps provided by th e 



WMAP scien c e team (IFosalba fc Szapudil (120041 ): iPatanchon et al.l ( 12005( 1: 



Eriksen et all (l2007al Jbl)). as well as using uncleaned maps where some mo dels 



of foregrounds and (or) detector noise are necessary (jEriksen et al.l (l2008blJah 1 . 



Other approaches for foreground cleaning, usi ng needlet coefficie nts ( iDelabrouille et al. 



( 120091 )) and harmonic variance minimization (IKim et al. (l2008al lblH have also 
been proposed. 

The organization of our paper is as follows. We describe the methodology 
of power spectrum estimation in Section [2j We present the power spectrum 
results along with their covariance structure in Section [31 Finally, we con- 
clude in Section [U 
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2 Method 



The basic procedure f or ext racti ng the temperatu re power spectrum is de- 
scribed in ISaha et al.l (120061 ) and ISaha et al.l (120081 ). Here we generalize this 



to include polarization. The basic maps for the case of polarization are avail- 
able in terms of the Stokes parameters Q and U. Since these are coordinate 
dependent quantities it is more convenient to w o rk wi t h the coordinate inde- 
pende n t E and B modes (jZaldarriaga &: Seljak! (119971 ) ; ISeljak fc Zaldarriaga 
(Il997h : iKamionkowski et all ril997t )).' Another problem with using Q and U 
maps i s that the E and B modes mix with one another when one applies a sky 



mask dJaffe et al.l (120001); iTegmark fc de Oliveira-Costal (120011 ): iBunn et al 



( 120031 ) ; ISmith fc Zaldarriagal (120071 )) to remove heavily contaminated Galac- 
tic regions. This demands an extra data processing step to isolate the actual 
CMB E and B mode power spectra from their mixture. To avoid this prob- 
lem we start by converting full sky Q and U maps to full sky E and B maps 
and apply mask when ever required on the resultant maps. This is similar 
to what is proposed in iBetoule et al.l (120091 ) for estimating r = T/S for the 
Planck satellite mission and the Experimental Probe of Inflationary Cosmol- 
ogy (EPIC). To obtain the E and B maps we first expand the full sky spin 
±2 fields (Q ± ill) in terms of spin-2 spherical harmonics ±2Yi m (n) 



(Q + iU)(n) = ^ a Vm 2 Y im{n) 

im 

(Q - iU){n) = ^2 a -2M -2Yi m (n) 



(1) 



im 



Since both Q and U are real, one can show that the expansion coefficients 
obey a*_ 2 i m = a 2,i-m- The spin-0 E and B are now obtained by the usual 
spherical harmonic transform, 



E{h) 
B(h) 



l>2,\m\<l 

Yl a lm Y lmin) 



(2) 



;>2,h<« 



where 



1 



Hm 



{&2,lm + a - 



21m) 



— {0>2,lm ~ Q>-2,lm) 



(3) 
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This gives us 10 different full sky maps for each of the E and B fields cor- 
responding to the 10 WMAP Differencing Assemblies (DAs). The 10 DAs 
are labeled as K, Ka, Ql, Q2, VI, V2, Wl, W2, W3, W4 corresponding to 
the five different frequency channels K, Ka, Q, V and W. We note that the 
bands Q, V and W have 2, 2, 4 DAs respectively. 

We first eliminate the hig hly contam i nated Galactic plane from all the 
10 DA maps using P06 mask ( Page et al. _( 2007h V This p rocedure is slightly 



different from that described in ISaha et al.l (120061 . 120081 ). In the latter case 



the authors cleaned the entire unmasked sky in nine iterations and also pro- 
duced a full sky foreground cleaned temperature maps. However to eliminate 
potential residua l foreground conta mination arising from the Galactic plane 
the KQ85 mask (IGold et al.l (120081 )) is applied before computing the power 
spectrum. In the present work, we apply the mask right at the beginning 
since we are interested in extracting only the power spectrum. To cross-check 
our one iteration method we also divide the EE maps in several parts and 
then perform foreground removal in the iterative approach. We find that, 
the final power spectrum of this method is similar to the one iteration case. 

We select different possible linear combinations of 4 maps out of the 
available DAs as described below. The entire set of linear combinations are 
listed in Table [H Each of these linear combinations independently lead to a 
clean map. 

The cleaning is accomplished independently for each I, by linearly com- 
bining these maps with weights, wf, such that the spherical harmonic com- 
ponents of the cleaned map are given by, 

y^r— • 



Clean 
hm 



0=1 



(4) 



Here n c is the total number of maps used for cleaning. In the present case of 
4 channel cleaning, n. 



, - 4. The fact or Bf is th e circu larized beam transform 

are 



function for the frequency band a ( jHill et al.l (120081 ) ). The weights 
chosen so as to minimize the total power subject to the constraint 



T 



W l e = e 1 W l =1, 
where eo is a column vector with unit elements 

/ 1 \ 



e 



(5) 



(6) 



V 1 / 
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and W; is the row vector 



This constraint is required so 



as to preserve the CMB signal. The weights are obtained using the em 



Teemark et al. ( 


2003 


); - 


?eemark & Efstathiou 


(1999)7 


Eriksen et al. 


Delabrouille & Cardoso! 


(2009^) 



2008); 



W, 



C, 



(7) 



We label the resulting cleaned maps as Ci and CAi where i = 1, 2, 24. 
Here the maps Ci use the DAs K along with possible combinations of DAs 
from the bands Q, V and W. Similarly the maps CAi include Ka instead of 
K. The entire nomenclature is listed in Table [TJ In the case of the W band 
we average over two DAs before we start the foreground cleaning. Hence in 
Tabled] the notation W12, for example, refers to the average of the DAs Wl 
and W2. This averaging is not essential to the procedure and one may also 
directly use the original WMAP DAs. However averaging leads to a reduced 
detector noise in each cleaned map. 

After obtaining the 48 cleaned maps we cross-correlate them in selected 
combinations in order to reduce the contribution due to detector noise. We 
cross-correlate all pairs of maps such that the two cleaned maps in each pair 
are formed by distinct DAs. This gives us 24 cross-correlated power spectra 
on the masked sky. All the possible cross-correlations are listed in Table [2J 

We convert each of the 24 masked sky power spectra into full sky esti- 
mates of the underlying CMB power spectrum using the mode-mode cou- 
pling m atrix corresponding to the P06 ma s k foll o wing the MASTER a p- 
proach flHivon et al) (120021 ): iHinshaw et all (l2003h : iTristram etaD t00$ \). 
We then remove beam and pixel effects from each of these 24 full sky power 
spectra. Our final EE power spectrum is simply an uniform average of these 
24 cross-spectra. We rely upon Monte Carlo simulations to compute the error 
bars as well as possible bias in the extracted power spectrum. 

The neighboring multipoles in the power spectrum become coupled since 
the spherical harmonics lose orthogonality on a masked sky. Hence, to obtain 
full information about the two-point correlation function of the resulting 
power spectrum one needs to construct the covariance matrix, 

(AC/ACV) = ((Cj - (G»(CV - (CV»> , 

We compute the covariance matrix by Monte Carlo simulations. The corre- 
lations can be minimized by suitably binning the power spectrum. We use 
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K + Ql + VI + W12 = d 
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K + Ql + V2 + W13 = d4 
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TV 1 A 1 1 "\ 7"0 1 WT\ A rii r 

K + Ql + V2 + W14 = do 
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Ka + (4I + V2 + W14 = LAlo 


1 1 ~\ 7"0 1 1UOO flic 

K + Ql + V2 + Wzo = do 


Ka + + V2 + W23 = UAlo 


K + Ql + V2 + W24 = C17 


Ka + Ql + V2 + W24 = CA17 


K + Ql + V2 + W34 = C18 


Ka + Ql + V2 + W34 = CA18 


K + Q2 + VI + W12 = C19 


Ka + Q2 + VI + W12 = CA19 


K + Q2 + VI + W13 = C20 


Ka + Q2 + VI + W13 = CA20 


K + Q2 + VI + W14 = C21 


Ka + Q2 + VI + W14 = CA21 


K + Q2 + VI + W23 = C22 


Ka + Q2 + VI + W23 = CA22 


K + Q2 + VI + W24 = C23 


Ka + Q2 + VI + W24 = CA23 


K + Q2 + VI + W34 = C24 


Ka + Q2 + VI + W34 = CA24 



Table 1: List of the different combination of the DA maps, used to obtain 
the final 48 cleaned maps, denoted by Ci and CA\ where i = 1, 2, . . . , 24. 



CI (g> CA12 


C2 ® CA11 


C3 g 


CA10 


C4 (g 


CA9 


C5 <g 


CA8 


C6 (8 


CA7 


C7 (g> CA6 


C8 <g) CA5 


C9 (g 


CA4 


C105 


D CA3 


Cll€ 


J CA2 


C12$ 


J CA1 


C13® CA24 


C14® CA23 


C15§ 


-) CA22 


C165 


) CA21 


C175 


D CA20 




5 CA19 


C19® CA18 


C20^ CA17 


C21§ 


5 CA16 


C22£ 


3 CA15 


C235 


D CA14 


C24$ 


3 CA13 



Table 2: List of all the 24 cross-power spectra using the 48 cleaned maps, 
CI, C2,..., C24 and CAl, CA2,..., CA24. 
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a binning identical to that used by the WMAP team. Let C a denote the 
binned power spectrum. Then the covariance matrix of the binned spectrum 
is obtained as 

(AC a AC a ,) = ((C a - (C a ))(C al - (C a ,))) . 

The standard deviation obtained from the diagonal elements of the binned 
covariance matrix gives the error-bars on the binned final spectrum. Since 
the cosmic variance of the CMB power spectrum decays as ~ 1/(2/ + 1), 
the diagonal terms in the above correlation matrix decay with increasing 
multipoles. For a visual comparison of correlation between different bins we 
define a correlation matrix, C aa ', of the binned power spectrum, where, 

= (ACqACgl) 

aa ' ^({AC a f){{AC a ,f)' 
All the elements of this matrix are bound to lie between [—1, 1]. 



3 Results 



3.1 Temperature Power Spectrum 

The temperature power spectrum for 5-year WMAP d ata is obtained using 
the same procedure as described in ISaha et al.l (120081 )). The entire sky is 
divided into 9 regions depending on the level of foreground contamination. 
The whole cleaning is done with the iterative method, starting from the 
dirtiest region. 

All the 24 cross power spectra for temperature anisotropy are shown in 
Fig. [TJ We form an uniform average power spectrum by averaging over these 
24 cross power spectra. While obtaining the cleaned power spectrum we use 
the KQ85 mask to remove the residual foreground contamination near the 
Galactic plane. 

Even after applying the KQ85 mask which also removes a circular region 
around each of the known point sources, residual unresolved point sources 
cause a significant contamination in this power spectrum. There have been 
several attem pts to measure unresolved poi n t sou rc e contamination in th e 



CMB maps ( iTegmark fc de Oliveira-Costal (11998h : iKomatsu et al.l (120081 ): 
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Figure 1: The 24, binned, TT cross power spectra obtained by using the 
WMAP 5 year data. All the different combinations of cleaned maps used are 
shown in the box. The average power spectrum along with error bars (blue 
points) is also shown. The red line joins the individual binned averages. 



Huffenberger et al.l (120061 )). We estimate unresolved point sour ce contami- 
nation in our final power spectrum using the m odel present e d bvlNolta et al. 
(120081 ) following an approach similar to that in ISaha et all (120061 . l2008h .~ 

We compute the bias in the extracted spectrum by performing 150 Monte 
Carlo simulations. First we generate synchrotron, free-free and thermal dust 
maps corresponding to different WMAP frequencies using the Planck Sky 
ModeQ (PSM). Although several options are available in the PSM to generate 
galactic emission (e.g. with or without spinning dust, with or without small 
scales added), the largest scales in temperature are strongly constrained by 
observation, and the impact of the choice of a particular model is not a major 
source of uncertainty. In our simulations, we use a single set of galactic 
emission maps, which comprise a two-component dust model based on SFD 
model 7, synchrotron map with varying spectral index in agreement with 
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|http://www.apc.univ~paris7.fr/ AP C _C S / Recherche /Adamis/PSM / psky-en . php 
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the first year WMAP data, and free-free emission with fixed spectral index, 
obtained from an H-alpha template corrected for galactic dust extinction. 
The exact polarisation properties of the galactic foregrounds, in particular 
that of dust emission, are poorly constrained by observations. For the present 
work, we use version 1.6.4 of the PSM (see Betoule et al. (2009) for details 
about the polarised galactic emission). In the next ste p we randomly genera te 
CMB maps assuming the standard LCDM model ( Spergel et al.l ( 2003 )). 
Each random realization of the CMB map is then added to the combined 
mixture of all three foreground components corresponding to the 5 WMAP 
frequencies. Using the 5 year beam transform functions for different DAs 
provided by the WMAP science team, we transform the 5 resulting maps 
into 10 maps. Each map at this step has a resolution appropriate for the 
corresponding DA. We then generate random noise maps corresponding to 
each detector. The random noise maps are generated by sampling a Gaussian 
distribution with unit variance and then multiplying each Gaussian v ariab le 
by Cq / \/Np, where do is the noise per observation ( IHinshaw et al. 
and N p the effective number of observations at each pixel. The values of <To 
depend on the DA, with the smallest value for the K band DA and largest 
for the W band DAs. Finally the noise maps are added to the CMB plus 
foreground maps for different DAs. These maps with CMB signal, detector 
noise and foreground are then passed through the same power spectrum 
estimation method as in the case of observed data. The mean of the 150 
extracted spectra gives the final simulated power spectrum. The standard 
deviation of the 150 simulations gives the error. The difference between the 
simulated power and the input LCDM power gives a measure of the bias in 
our method. This bias is subtracted from the extracted power spectrum in 
order to get the final result. 

The precise magnitude of the bias depends on the theoretical model with 
which we compare our extracted power spectrum. In other words, before 
we compare our extracted power to a theoretical model, we must correct 
for bias using the corresponding model power spectrum. Here we use the 
WMAP best fit LCDM model to compute the theoretical power spectrum. 

The final temperature power spectrum using the Internal Power Spectrum 
Estimation (IPSE) method, after correcting for bias, is shown in Fig. [2j We 
find that it is in good agreement with the WMAP best fit LCDM model. The 
simulation results are also shown in Fig. [2j After bias correction we find the 
quadrupole power (1(1 + l)Cj/27r) equal to 532 fiK 2 compared to the value of 
213.4 fj,K 2 estimated by the WMAP science team. The quadrupole extracted 
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from IPSE is, therefore, in much better agreement with the theoretical model. 

The correlation matrix, Eq. [SJ for temperature power spectrum is shown 
in Fig. [3j We find, as expected, that the off diagonal matrix elements are 
negligible compared to the diagonal elements. 
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Figure 2: The final, binned, TT power spectrum using IPSE (red dots) along 
with error bars for the 5 year WMAP data, after subtracting the bias ex- 
tracted using simulations. The error bars are also obtained from simulations. 
The WMAP LCDM best fit model (black line) and the simulation results 
(green crosses) are shown for comparison. The bottom panel shows the cor- 
rection made for residual power from unresolved point source contamination. 



3.2 The TE power spectrum 

The WMAP polarization CMB maps are cleaned using a single iteration 
rather than the nine iteration procedure followed for temperature anisotropy. 
We apply the P06 mask right in the beginning. The mask removes 27% of the 
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Figure 3: The correlation matrix elements, C aa > } defined in Eq. [SJ for the 
temperature power spectrum, plotted with respect to the bin indices a, a'. 



entire sky region near the galactic plane. The cleaning algorithm is applied 
only to regions outside the P06 mask. Hence we make no attempt to produce 
a full sky cleaned polarization map. 

The error bars plotted in the TE power spectrum are obtained by Monte 
Carlo simulatio ns. We gene r ate 1 50 random samples of data using the 
WDM model flSpergel et all (l2003f n and simulated foregrounds and detec- 



tor noise. First we generate the synchrotron and thermal dust polarized 
foreground maps corresponding to different frequencies in terms of Q and U 
maps using the PSM version 1.6.4. The free-free emission is not polarized 
and hence not included. The anomalous dust emission is also assumed to be 
unpolarized and thus is not included either. Using the HEALPix[f| command 
synf ast we generate random realization of CMB polarization maps in terms 
of Q and U maps. The random CMB realization and foreground maps are 
smoothed by the beam functions corresponding to ten different DAs. We 



2 http: / /healpix.jpl. nasa.gov 
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next obtain E-mode polarization maps from these Q and U maps. Then we 
generate random noise maps for each DA in terms of Stokes parameters Q and 
U using Cholesky decomposition technique for generating correlated Gaus- 
sian random variables using the WMAP supplied 2x2 QU intra-pixel noise 
covariance matrices. These Q and U noise maps are converted to E mode 
noise maps. The final E-maps including detector noise, foregrounds and 
CMB signal are passed through the same cleaning pipeline as the observed 
polarization data. In order to minimize the correlation among neighboring / 
modes, the final power spectrum is binned in the same way as the WMAP 
5 year result. Here also the standard deviation obtained from the diagonal 
elements of the binned covariance matrix is used as the error bars on the 
binned final spectrum extracted from the WMAP data. 

The extracted TE power spectrum along with the WMAP results and the 
best fit LCDM model is shown in Fig. HI The binned TE power spectrum, 
using the same binning scheme as used by the WMAP team, is shown in Fig. 
|5j The error bars are computed by simulations. The simulation results are 
shown in Fig. EJ We find that the bias is small for all the bins. Only at 
small I, I < 10, do we find a noticeable negative bias. For larger I, the bias is 
practically negligible. The bias corrected TE power spectrum is shown in Fig. 
[3 The spectrum obtained by the WMAP science team as well as their best 
fit LCDM model is also shown. We find good agreement with the WMAP 
result. However we obtain slightly smaller error bars. We discuss the possible 
cause for this in the next subsection. The correlation matrix elements, Eq. 
IH1 are shown in Fig. [HJ Here also we see that they are dominated by diagonal 
elements. 

3.3 The EE power spectrum 

The binned EE power spectrum, using the WMAP binning procedure, along 
with the simulation results, is shown in Fig. [9] for the 5 year data. In this 
figure we also show the EE power spectrum extracted by the WMAP team 
along with the WMAP best fit LCDM model. In Fig. |we follow exactly 
the WMAP binning procedure and consider only the multipoles Z > 50. The 
results for low multipoles / < 50 are shown separately. We find that our 
extracted power spectrum is in good agreement with that obtained by the 
WMAP team but with significantly smaller error bars. Furthermore the 
binned simulated power spectrum is found to be close to the input LCDM 
power over the entire multipole range. Only at small I do we find a significant 
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Figure 4: The cleaned TE power spectrum using IPSE (solid blue line) along 
with WMAP result (dashed black line). The WMAP best fit LCDM power 
spectrum (thick solid red line) is shown for comparison. 

positive bias. For the remaining multipoles the simulation results match 
the input power within error bars. The bias corrected power spectrum is 
shown in Fig. [TU1 We see that the bias corrected spectrum is in reasonable 
agreement with the best fit LCDM model. The correlation matrix elements 
are shown in Fig. [HJ We again find that the correlation matrix is dominated 
by diagonal matrix elements. 

Our extracted power spectrum along with the simulation results for low 
I (I < 50) are shown in Fig. [T5J The WMAP power spectrum as well as 
their best fit LCDM model is also shown. Here we have chosen the binning 
that was used by the WMAP team for the TE power spectrum. The bias 
corrected power spectrum is shown in Fig. [131 

Since the power spectrum estimation methodology described in this work 
is fundamentally different from WMAP team's approach we should not ex- 
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Figure 5: The final binned TE power spectrum using IPSE (blue diamonds) 
along with error bars, compared with the WMAP results (red dots). The 
WMAP best fit LCDM result (solid pink line) is also shown. 

pect both methods to produce identical error-bars on the derived power spec- 
tra, although the power spectra themselves, obtained by using the two meth- 
ods, are in reasonable agreement with each other. We find that our error-bars 
on the polarization power spectra are smaller compared to those obtained by 
the WMAP science team. This effect could be explained by noting that 
we use more detector maps for polarization power spectrum estimation than 
WMAP science team. Using more detectors increases the signal to noise ratio 
of the cleaned map by decreasing the effective noise level. The reduced noise 
level leads to lower error-bar on our polarization power spectra. Further- 
more, in the case of noisy polarization data the weights tend to be inversely 
proportional to the noise. Due to this inverse noise weighting each cleaned 
map has lesser noise in comparison to the least noisy K or Ka band maps. 
This may be another reason for smaller error bars in the case of polarization 
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Figure 6: The binned TE power spectrum, along with error bars, obtained 
from simulations (red dots). The input LCDM model (solid black line) is 
also shown. 



power spectrum. 



4 Conclusion 

We have used a model independent method to estimate the CMB temperature 
and polarization power spectrum using WMAP 5-year data. The method is 
based on the assumption that the CMB signal is independent of frequency 
in thermodynamic temperature unit. Since the foregrounds are frequency 
dependent in this unit, it is possible to minimize the foreground power by 
making a linear combination of CMB maps with a suitable choice of weights. 
For foreground minimization we use the CMB maps in harmonic space. For 
the case of polarization, the raw full-sky Q and U maps are first converted 
to E and B maps to avoid any mixing of E and B modes. The total number 
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Figure 7: The bias corrected, binned TE power spectrum using IPSE (blue 
diamonds) along with error bars, compared with the WMAP results (red 
dots). The WMAP best fit LCDM power spectrum (solid pink line) is also 
shown. 

of maps available for each field for WMAP data is ten, corresponding to the 
ten DAs. We create several cleaned maps by choosing different subsets of ten 
DAs, such that each set contains only four DAs. The detector noise power 
is minimized by cross correlating cleaned maps obtained from distinct DAs. 
This leads to a considerable reduction of the detector noise power since, to 
a good approximation, noise is uncorrelated among different detectors. 

By utilizing all the ten WMAP DA maps to estimate the polarization 
power spectra, we are able to provide more stringent constraints on the spec- 
tra in comparison to that obtained by the WMAP science team. We find 
that the error-bars on our TE and EE power spectra are smaller than those 
obtained by the WMAP science team on all angular scales. Another possible 
reason for why we get smaller error bars in the case of polarization power 
spectra is that in a noisy data the weights tend to combine the maps in in- 
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Figure 8: The covariance matrix elements, C aa / } defined in Eq. [5J for the 
TE power spectrum plotted with respect to the bin indices a, a'. 



verse noise weighted manner. This results in each cleaned map having lesser 
noise than the least noisy K or Ka band maps. 

In the case of TT power spectrum we find that our procedure does not 
remove all the unresolved point source contamination. This contamination 
is significant at small angular scales where the detector noise is also very 
large. Hence here our internal cleaning is not very efficient. The residual 
unresolved point source contamination is removed b y using the WMAP point 
source model, as described in detail in ISaha et al.l (1200a [20081) - 

We have performed detailed simulations of the TT, TE, EE power spec- 
trum, using the WMAP best fit LCDM model, along with foreground and 
detector noise models, in order to determine if there exists any bias in the 
extracted power. In all cases the bias is found to be small for the entire mul- 
tipole range. The extracted power, with or without bi as correction, is fo und 
to be in good agreement with the WMAP results. In ISaha et al.l (120081 ) the 
authors noticed a negative bias at low / in the TT power spectrum. The neg- 
ative bias arises due to a chance correlation between the CMB and the fore- 
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Figure 9: The CMB polarization EE power spectrum using IPSE (blue dia- 
monds) along with error bars, compared with the WMAP result (red dots). 
The WMAP best fit LCDM power spectrum (pink solid line) and the sim- 
ulation results (black crosses) are also shown. 



grounds. After correcting for the negative bias, we find that the quadrupole 
for the WMAP 5 year data shows much better agreement with the LCDM 
model, in comparison to the result obtained by the WMAP science team. 
Excluding / = 2, we find negligible bias at all multipoles except at very large 
I values, where we find a small positive bias. For the case of TE power spec- 
trum we also find a small negative bias at low I, I < 10. For larger I values 
the bias is negligible. For EE power spectrum also the bias is found to be 
small compared to the corresponding error bars. A significant positive bias 
is found only at low /. 

To summarize, we have performed a completely independent reanalysis 
of WMAP 5 years temperature and polarization data. Our procedure uses 
primarily the CMB data. Hence it is free from any bias that might result 
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Figure 10: The CMB polarization EE power spectrum using IPSE (blue 
diamonds) along with error bars, after correcting for bias. The WMAP 
result (red dots) and the theoretical LCDM spectrum (pink solid line) are 
shown for comparison. 



from the inadequacies and inaccuracies of the foreground modeling. The 
foreground templates and detector noise modeling is utilized only for the 
purpose of bias analysis and error estimation. The bias is found to be small 
for all the spectra over the entire multipole range. Our results verify the basic 
power spectra results obtained by the WMAP Science team. We expect that 
the method will be very useful for analyzing data from future CMB probes 
such as Planck. 
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Figure 11: The correlation matrix elements, C a>a i, defined in Eq. El for the 
EE power spectrum plotted with respect to the bin indices a, a'. 
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